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We consider geometrical clusters (i.e. domains of parallel spins) in the square lattice random 
field Ising model by varying the strength of the Gaussian random field, A. In agreement with the 
conclusion of previous investigation (Phys. Rev. E63, 066109 (2001)), the geometrical correlation 
length, i.e. the average size of the clusters, £, is finite for A > A c ~ 1.65 and divergent for A < A e . 
The scaling function of the distribution of the mass of the clusters as well as the geometrical 
correlation function are found to involve the scaling exponents of critical percolation. On the other 
hand the divergence of the correlation length, £(A) ~ (A — A c ) _ ", with v m 2. is related to that 
of tricritical percolation. It is verified numerically that critical geometrical correlations transform 
conformally. 



I. INTRODUCTION 

The random-field Ising model (RFIM) is a prototype 
of random systems in which the disorder is coupled to 
the order-parameter of the system*. It has experimen- 
tal realizations, such as a diluted antifcrromagnct in a 
fielc£. For some time there was a debate about the lower 
critical dimension, di, of the systems. While domain- wall 
stability arguments by Imry and Ma ? predicted di — 2, 
perturbative field-theoretical calculations^ have lead to a 
different value of di = 3. Later, Fisher argued^ that due 
to the existence of several metastabile states the per- 
turbative renormalization does not work. Indeed, exact 
results by Bricmont and Kupiainen- show that in the 3d 
RFIM there is ferromagnetic order and later Aizenman 
and Wehr— has proven rigorously that at d — 2 the Gibbs 
state is unique, thus di = 2. 

Although there is no ferromagnetic long-range order in 
the RFIM at d = 2 there are several interesting questions 
which are related to the structure of clusters of parallel 
spins and the corresponding geometrical correlations in 
the system. Even at zero temperature, at T = 0, the 
ground state of the system is not trivial. This is the re- 
sult of a competition between the ordering effect of the 
ferromagnetic nearest-neighbor interaction, J, and the 
disordering effect of the random field with a strength, A, 
which is given by the variance of the distribution. In the 
limit of strong random fields, A/J^S> 1, the direction of 
the spins follows the actual direction of the random fields 
and the ground-state is equivalent to a site-percolation 
problem^ in which an up (down) spin corresponds to an 
occupied (empty) lattice site. Since the occupation prob- 
ability, p = 0.5, is below the site-percolation threshold 8 , 
p c = 0.593, in the ground state the domains of the par- 
allel spins have only a finite extent, the linear size of 
which, £, is given by the correlation length of percolation 
at p = 0.5. As the strengths of the random field is de- 
creased there is a tendency of the formation of larger par- 
allel domains, the typical size of which is a monotonously 
increasing with decreasing A. It is an interesting ques- 
tion, if £ stays finite for any A > 0, or there is a finite 



limiting value, A c > 0, at which £ becomes divergent. 
According to recent numerical workS this second scenario 
holds, so that for weak enough random fields, A < A c , 
the clusters of parallel spins, which are called geometrical 
clusters, percolate the sample. 

Geometrical clusters in non-random spin systems, such 
as in the Ising and the Potts models, have been defined 
for a long timeiS and their properties have been inten- 
sively studiedii at a finite temperature, T > 0. In two 
dimensions the geometrical clusters percolate at a tem- 
perature, T g , which corresponds to the critical point of 
the system^: T g = T c . In 2d their fractal dimension 
can be obtained through conformal invariancei^ and this 
value is generally different from the fractal dimension of 
the so called Fortuin-Kasteleyn clusters^*. The Fortuin- 
Kasteleyn clusters are represented by graphs of the high- 
temperature series expansion of the models and the cor- 
responding fractal dimension is directly related to the 
scaling dimension of the order parameter. In three di- 
mensions the geometrical clusters of any spin orientations 
percolate in the complete paramagnetic phase. Here the 
percolation transition temperature, T g , is defined for the 
minority spin orientation, so that for T < T g < T c the 
geometrical clusters of minority spins do not percolate^. 
In this respect geometrical clusters play a somewhat anal- 
ogous role in the 2d RFIM at T — (by varying A) as 
the geometrical clusters of minority spins in the 3d Ising 
model for T < T c (by varying T). 

In this paper we consider the properties of the geo- 
metrical clusters in the 2d RFIM. In previous numerical 
works a homogeneous field, H, was also applied and the 
behavior of the magnetization and the susceptibility is 
studied by varying H and A. The obtained schematic 
phase-diagram of the model is presented in Fig^ which 
will be discussed in details in Sec[n] In another worki& 
nonequilibrium critical relaxation of the 2d Ising model 
has been studied, in which the initial state was prepared 
as the ground state of the RFIM. By varying A and set- 
ting H = the geometrical phase transition is found to 
influence the properties of the nonequilibrium dynami- 
cal processes by introducing new finite time and length 
scales. 
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In the present paper we study those aspects of the ge- 
ometrical clusters in the RFIM which have not been yet 
investigated in previous work. Here we set H = 0, vary 
the strength of disorder and study the distribution of the 
mass, M of the geometrical clusters, R(M, L) in a finite 
system of linear size, L. In particular from the scaling 
form of R(M, L) we calculate both the fractal dimension, 
df, and the distribution exponent^., r. We also measure 
the geometrical correlations, G(r). This is the average 
over all spin pairs of distance, r, with a contribution one, 
if both spins are in the same geometrical cluster and zero 
otherwise. We show that G(r) plays the role of a stan- 
dard correlation function in the vicinity of a second-order 
phase transition point. The correlation length associated 
to G(r), £(A), is found to be finite for A > A c and di- 
vergent for A < A c , thus defining two different regions, 
see in FigQ] In the percolating regime, A < A c , G(r) 
exhibits quasi-long-range order, G(r) ~ r~ v . We have 
also studied conformal aspects of G r in the percolating 
regime. 

The structure of the paper is the following. The model, 
the numerical method and the known results are given in 
Sec[n] Our results are presented in Sec lIIII and discussed 
in SeclTvl 



II. MODEL, NUMERICAL METHOD AND 
KNOWN RESULTS 



The RFIM is defined by the Hamiltonian: 



(1) 



where <7j = ±1 is an Ising spin located at site i of a square 
lattice. The ferromagnetic nearest-neighbor coupling is 
chosen to be J = 1, whereas the hi random fields are 
independent variables which are taken from a symmetric 
Gaussian distribution: 



P(hi) = 



1 
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which has a variance A. For completeness in Eq.JSJ we 
have also included a homogeneous field of strength, H. 
During our numerical calculations we shall put H = 0. 

We are interested in the properties of the system at 
T = when all information is encoded in its ground 
state. The problem of finding the ground state of the 
RFIM is a non-trivial task, but can be exactly solved 
through a mapping to a maximum flow problemit For 
this latter problem there are very efficient combinatorial 
optimization algorithms 1 ^ which work in strongly poly- 
nomial time. With this method finite systems as large as 
L = 500 can be treated numerically. 

At H = and for A > in the ground state of the 
RFIM there is no magnetic long-range orderi This can 
be illustrated by the average spin-spin correlation func- 
tion, C(r) — (ao<7 r ), which is averaged over all pairs of 



spins having a distance r apart and (...) stands for the 
calculation in the ground state. C(r) goes to zero expo- 
nentially with the distance, r. In the presence of a ho- 
mogeneous field, H ^ 0, there is a finite magnetization, 
m(H) = (a) =/= 0, and one can measure the magnetic sus- 
ceptibility, x{H)- By varying however, x{H) displays 
no singularity^. 

The ground state of the RFIM can be visualized by 
denoting the up and down spins by different symbols. In 
this picture the domains of parallel spins form geometri- 
cal clusters, the structure of which depends on the value 
of A (and also H). The typical size of the clusters (of 
majority spins), £, serves as a bulk length-scale in the 
problem. Between large, oppositely magnetized clusters 
there is an interface, which is smooth for small lengths, 
up to ( < lb, and its roughness is seen only for larger 
scales, for I > If,. Here 4 is the so called breaking-up 
length. For weak random fields the breaking-up length 
asymptotically behaves asi^ 



l b ~ exp(A/A 2 



(3) 



thus it is divergent as A — > 0. The existence of the 
breaking-up length imposes limitations on the possible 
numerical calculations. The linear size of the system, L, 
should be much larger than lb, therefore one can not use 
a too small value of A. Due to this reason we have not 
reduced A below 1.2 — 1.4 for the system sizes we could 
treat numerically. 
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FIG. 1: (Color online) Schematic phase-diagram of geomet- 
rical clusters of the 2d RFIM as obtained in Ref . 9] . In the 
regions, F f and F j, the majority spin clusters are non- 
percolating, in P \ and P j, they are percolating. The phase 
boundary is indicated by dashed lines. The transition along 
arrow (1) is in the standard percolation universality class. 
Investigations in this paper are restricted to H — 0, in par- 
ticular in the vicinity of the multicritical point along arrow 
(2). 

Some aspects of the structure of the geometrical clus- 
ters for a non-zero homogeneous field has been studied 
in Ref.jjJ, see the phase-diagram in Fig^ For a fixed 
A > A c the size of the majority spin clusters, £, is finite 
for small H , but monotonously increasing with increas- 
ing H. There is a threshold value, H p (A), when the 
majority spins percolate the sample and for H > H p (A) 



a finite fraction of the majority spins belongs to the infi- 
nite cluster. Close to the percolation point, see arrow (1) 
in FigQ] the percolation probability is found to depend 
on the scaling combination: L(H p — H)"*, with v p = 4/3 
characteristic of the 2d short range percolation transition 
point 8 . Also the fractal dimension of the geometrical 
clusters at the percolation point is found^ to be com- 
patible with that of standard percolation^, df = 91/48. 
According to the numerical results as A is decreased H p 
tends to zero at A c , as H p ~ (A — A c )*, with A c w 1.65 
and (j) « 2.05. 

In the following section we concentrate on the proper- 
ties of the geometrical clusters at H = 0. In particular 
we study the distribution of their mass and introduce and 
investigate a correlation function, which is associated to 
geometrical clusters. In these investigations we approach 
the multicritical point at A = A c and H = along arrow 
(2) in Fig^ With these investigations we want to shied 
light to the possible geometrical phase transition of the 
system which takes place at A c . 



III. RESULTS 

We have considered the RFIM on L x L square lat- 
tices with periodic boundary conditions in both direc- 
tions. The strength of the random field is varied in the 
range 1.2 < A < 4., where the lower limiting value is 
set by the condition, L ^> Zh(A), where the breaking- 
up length is given in Eq.(PJ. Using combinatorial opti- 
mization algorithm we have calculated the ground state 
configuration exactly for different finite samples up to a 
linear size L — 500. Averages are performed over typi- 
cally 10000 different realizations of disorder. 



A. Distribution of the mass of geometrical clusters 

We start to study the cumulative distribution of the 
mass (number of spins), M, of the clusters, which is de- 
noted by R(M, L). This measures the fraction of clusters 
having at least a mass of M. According to scaling the- 
ory developed for standard percolation^ this distribution 
asymptotically behaves as: 

R(M, L) = Ar T R{M/L df ) (4) 

with a characteristic exponent, r = 2/df. Thus the only 
parameter which enters into this expression is the fractal 
dimension of the clusters, df. This relation is expected to 
hold at the critical point, where the correlation length, £, 
associated to geometrical correlations is divergent. How- 
ever, Eq.Q can be a good approximation outside the 
critical point, too, provided the correlation length, £, is 
much larger than the size of the system, L. To compute 
R(M, L) we have performed cluster statistics over 10000 
different random samples and the distributions obtained 
for different finite systems are scaled together using the 
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FIG. 2: (Color online) Scaling plot of the cumulative distribu- 
tion of the mass of the geometrical clusters at A = 1.8, using 
the fractal dimension of standard percolation. In the inset 
a) the effective fractal dimension is presented as calculated 
from the optimal scaling collapse of the curves for different 
values of A. In the inset b) extrapolation of the r exponents 
calculated at different finite systems is presented for A = 1.4 
(x), A = 1.6 (+) and A = 1.8 (*). 



relation in Eq.QJ. As an illustration in Fig^Jwe present 
the scaled cumulative distribution functions at the spe- 
cific value, A = 1.8, in which we used the fractal dimen- 
sion of standard percolation 8 , df — 91/48. The scaling 
collapse in Fig|2 is indeed satisfactory. Repeating the 
same procedure in the range of 1.4 < A < 2.6 we have 
obtained still acceptable data collapse. This means that 
in this regime the condition £ L is satisfied and we 
can conclude that the critical point, A c , is also located 
in this range. 

We have also estimated the value of the fractal dimen- 
sion from the optimal scaling collapse of the curves. For 
this we have measured the surface of the overlap of the 
scaled curves for different values of the parameter, df, 
and at the true fractal dimension the overlap surface is 
minimal, see in Refi2&. The measured effective, i.e. A- 
dependent fractal dimensions are presented in the inset 
a) of FigEl In the range of 1.6 < A < 2. the effec- 
tive fractal dimensions are approximately constant and 
consistent with df = 1.89(2). For A > 2. the effective 
df start to decrease, which is a clear indication that the 
correlation length is not large enough. The decrease of 
df for A < 1.6, is probably due to the large value of the 
breaking-up lengths, which starts to approach the size of 
the system. 

From the slope of the curves in Fig[3in the log-log plot 
we have measured the exponent, r, which is defined in 
Eq.|@J. The measured effective exponents for different 
finite sizes are plotted in the inset b) of FigEl f° r three 
different values of A. Here strong finite-size dependence 
is visible so that we have to perform an extrapolation 
procedure in terms of 1/L. The extrapolated value, as 
seen in the inset b) of Fig[3 is given by r = 1.055(3), 
and this is independent of A we considered. We note 
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that our estimate fits very well to the theoretical value of 
standard percolation^: t = 96/9f = 1.055. We can thus 
conclude that the distribution function of the mass of the 
geometrical clusters nicely satisfies the scaling prediction 
in Eq.(£lJ and the fractal dimension coincides with that 
of standard percolation. 

B. Geometrical correlations 

Here - using the analogy with ordinary percolation - 
we introduce the geometrical correlation function, G(r), 
which is associated to the geometrical clusters. By def- 
inition between two spins the geometrical correlation is 
one, if both spins belong to the same geometrical cluster 
and zero otherwise. G(r) is then obtained by averaging 
the correlations over all pairs of spins having a distance, 
r. In the non-critical state, where £ < oo, the geomet- 
rical correlation function is short ranged and expected 
to decay asymptotically as G(r) ~ exp(— r/£). On the 
other hand at the critical point the geometrical correla- 
tion function has an algebraic decay, G(r) ~ r~ v , where 
the decay exponent, t], is related to the fractal dimension 
as, r\ = 2(d — df), through scaling theory. 

In practical calculations the system has a finite size, L, 
and in the critical region the geometrical correlation func- 
tion, G{r,L), is expected to behave asymptotically as: 
G(r,L) = r~ n G(r/L). Here the scaling function, G(y), 
should approach a limiting finite value for small argu- 
ments, so that the decay exponent, 77, can be calculated 
for large finite systems, if 1 <C r <C L. As an illustration 
we shown in FigE] the geometrical correlation function 
for different finite systems in log- log plot for A = 1.8. 
Indeed the slope of the curves (of the largest systems) in 
the region 10 < r < 40 is compatible with the theoretical 
prediction of ordinary percolation: r] = 5/24 = 0.208. 



0.5 - 




10 



FIG. 3: (Color online) Geometrical correlations in a log-log 
plot for A = 1.8 and for different system sizes. The theoret- 
ical prediction of the asymptotic behavior is indicated by a 
straight dotted line. 

In order to obtain a more accurate estimate of the 



asymptotic behavior of the geometrical correlation func- 
tion we measured the correlations between two points 
of maximal distance, i.e. at r = L/2. The maxi- 
mal distance correlation function defined in this way, 
G(L) = G(r — L/2,L), is expected to behave in the 
critical region as: G(L) = L _r, G(L/^), thus it is advis- 
able to consider the scaled maximal distance correlation 
function, L n G{L) = G(L/£). In FigHwe have plotted 
this function for different values of A. 
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FIG. 4: (Color online) Scaled maximal distance geometrical 
correlations for different values of A. The scaled correlations 
seem to reach a constant limiting value for A < 1.6, thus here 
the correlation length is divergent. For stronger random fields 
the correlations decay to zero, although the correlation length 
is still a large finite value. Inset: the correlation length as a 
function of A — A c in a log-log plot. The best fit is obtained 
with A c = 1.65 and with an exponent, v = 1.98(5). 

Here one can identify two different regimes as far as 
the large L dependence of the scaled correlation func- 
tion is concerned. For strong random fields, A > A c , the 
scaled correlations start to decay to zero, so that the cor- 
relation length is finite and the geometrical clusters are 
non-percolating. On the other hand for weak random 
fields, A < A c , the scaled correlations approach a finite 
limiting value, so that the correlation length is divergent 
and we are in the phase in which the geometrical clusters 
are percolating. This conclusion is the same as that ob- 
tained through the analysis of the spanning probability 
in Ref A 

The threshold value of A c can be calculated by 
analysing the behavior of the correlation length, which 
can be deduced from the asymptotic dependence: 
G(L/£) ~ exp[— L/(2£)]. Close to the critical point 
the correlation length is expected to be in the form: 
£(A) ~ (A - A c )-". As illustrated in the inset of Figgl 
the best fit is obtained with A c 1.65, as in Rcf. 9 , 
whereas for the correlation length critical exponent we 
have obtained: v — 1.98(5). This exponent describes the 
singularity of £ along arrow (2) in Fig|T] If we want to 
calculate the divergence of the correlation length at A c 
but with a small homogeneous field, H, i.e. along arrow 
(3) in Fig^ we make use of the known result about the 
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critical boundary, H p ~ (A — A c )^. Than from scal- 
ing theory we obtain £(H p ,A c ) <~ £(0, A) ~ H^", with 
v = u/cj) = 0.97(5). We note that this value corresponds 
to the inverse of the second thermal exponent of tricriti- 
cal percolation^: l/yt2 = 1- 

C. Geometrical correlations in the strip geometry 




1/L w 

FIG. 5: (Color online) Effective decay exponents calculated 
in periodic strips of width L w through Eq.Q. The dashed 
straight line indicates extrapolation through 1/L W . Inset a): 
Geometrical correlation function at A = 1.8 in a strip of width 
L w = 40 and for different lengths, L\. In a semi- log plot the 
slope of the limiting curve is indicated by a dashed line and 
given by l/£(L m ). Inset b): A-dependence of the effective 
decay exponent for L w = 40. 

Critical correlations are generally invariant under 
conformal transformations and conformal invariance 
has very important consequences in two dimensional 
systems*^ One important application of the conformal 
method is to transform the correlation functions from one 
geometry to another. In this respect it is of interest if also 
geometrical correlations are conformally invariant and if 
one can transform them between different geometries. To 
answer to this question we have studied geometrical cor- 
relations in the RFIM also in the strip geometry. It is 
known that the infinite plane, described by the complex 
variable, z — x + iy, can be mapped into a periodic strip 
of width, L Wl having the complex variable, w = u + iv, 
through the logarithmic conformal transformation: 

w = — — In z . (5) 

Then critical correlations in the plane, G(z) ~ \z\~ v , are 
transformed into an asymptotic exponential decay along 
the strip: 

G Li » ~ expHV£(^)) » (6) 

where the correlation length, £(L W ), is proportional to 
the width of the strip and given by: 

Z(L W ) = ^ . (7) 
irr] 



Thus the proportionality factor contains the decay expo- 
nent in the plane, therefore Eq.J7|) is called the correla- 
tion length - exponent relation^S. 

To check the validity of the correlation length - ex- 
ponent relation for the RFIM we used strips of widths, 
L w = 16, 20, . . . , 40 and calculated the geometrical corre- 
lation function along the strip. The lengths of the strips 
are taken so large, that the calculated exponential decay 
in Eq.© becomes independent of it. This is illustrated 
in the inset a) of Fig[SJ Generally we went at least up to 
a length of 16 x L w sites. After measuring the correlation 
length for a given L w we have calculated from Eq.Q) an 
effective, i.e. L^-dependent exponent. This calculation 
is then repeated for several values of A in the range of 
1.4 < A < 2.0. For a given L w the A-dependence of 
the exponents are found to be smaller than the actual 
error of the calculation, as shown in the inset b) of Fig|5] 
Therefore in Fig|S] we have plotted an average, i.e. A- 
independent exponent as a function of l/L w . As seen 
in this figure the effective exponents have a size depen- 
dence and an extrapolation through 1/L W leads to the 
estimate: rj = 0.210(4). This value is in good agreement 
with the conformal prediction: ?/ = 5/24 = 0.208. Thus 
we conclude that our numerical study has confirmed that 
the correlation length - exponent relation is valid for ge- 
ometrical correlations in the RFIM. 



IV. DISCUSSION 

In this paper we have considered the RFIM in the 
square lattice and studied the structure of its ground 
state at T = 0. Since in the ground state of the sys- 
tem there is no ferromagnetic long-range order we have 
focused on the properties of the geometrical clusters con- 
sisting of parallel spins. In particular we have investi- 
gated the distribution of the mass of the geometrical clus- 
ters and studied geometrical correlations, which are as- 
sociated to geometrical clusters. In agreement with pre- 
vious investigations- the geometrical clusters are found 
to have a finite extent for strong random fields, A > A c , 
and being percolating for weak random fields, A < A c . 
For A < A c the fractal dimension of geometrical clus- 
ters is calculated by different methods and its value is 
found to be in good agreement with short range percola- 
tion. The geometrical correlation function of the system 
is calculated both in the plane and in the strip geome- 
tries. It is shown that the critical geometrical correlation 
function satisfies both scaling and conformal invariance. 
In the vicinity of the geometrical critical point the corre- 
lation length is found to diverge as: £(A) ~ (A — A c )~" , 
with an exponent, v w 2, whereas at the critical disor- 
der, A = A c , the divergence of £(H) in the presence of 
a homogeneous field involves the exponent, v ss 1. This 
latter is characteristic for tricritical percolation. 

We note that our study of geometrical clusters in the 
RFIM is somewhat related to the properties of the ran- 
dom bond Potts model (RBPM) in the large-g limit. As 
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was shown by Cardy and Jacobsen24 the interface Hamil- 
tonian of the RFIM separating different types of geo- 
metrical clusters, and the interface Hamiltonian of the 
RBPM separating different types of Fortuin-Kasteleyn 
clusters, can be mapped into each other. In this respect 
the absence of long-range order in the RFIM corresponds 
to the absence of phase coexistence in the RBPM, thus 
no first-order phase transition in the presence of bond 
disorder in 2d 7 . Although the interface Hamiltonian of 
the two problems are equivalent the cluster structure, 
in particular the value of the fractal dimension is differ- 
ent. For the RBPM the conjectured valuo 25 i 26 for the 
Fortuin-Kasteleyn clusters, dj = (5 + \/5)/4 = 1.809, is 
considerably smaller than that of the RFIM. Thus argu- 
ments based on the interface Hamiltonian can be used to 



study the stability of some phases, but the actual cluster 
structure is governed by such critical fluctuations which 
are encoded in the details of the Hamiltonian. 
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